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Tj" . I. INTRODUCTION 

o 

In the inspired final chapter of their classic book Misner, Thorne and Wheeler state that there are three levels 
of gravitational collapse. The first two of them are the gravitational collapse of the whole universe during the final 
stages of its recontraction, and the gravitational collapse of a star when a black hole is formed. The third level of 
gravitational collapse is the quantum fluctuation of spacetime geometry at the Planck scale of distances. To rephrase 
Misner, Thorne and Wheeler, "collapse at the Planck scale of distances is taking place everywhere and all the time in 
Q\ • quantum fluctuations in thegeometry and, one believes, the topology of space. In this sense, collapse is continually 
t-H j being done and undone..." 

The picture given by Misner, Thorne and Wheeler of the Planck scale physics in these sentences is very charming. 
It immediately brings into one's mind mental images of tiny wormholes and black holes furiously bubbling as a sort 
. of spacetime foam. This is indeed a wonderful picture but unfortunatelly the ideas it suggests have never been taken 
^vq ' very far. Instead of developing models of spacetime where spacetime consists of tiny wormholes and black holes the 
reseachers in the field of quantum gravity have gone into another direction and developed, among other things, a very 
f — • succesfull theory of loop quantum gravity which treats spacetime as a spin network constituted of tiny loops. 

The purpose of this paper is, in a certain sense, to revive some aspects of the old picture of spacetime as a foam 
of tiny wormholes and black holes. Actually, there are reasons to believe that at the Planck length scale microscopic 
black holes might indeed play a role in the structure of spacetime. For instance, suppose that we want to localize a 
Q . particle within a cube with an edge length equal to one Planck length 
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& Z P1 :=^«1.6-10- 35 m. (1.1) 

It follows from Hcisenberg's uncertainty principle that in this case the momentum of the particle has an uncertainty 
Ap ~ h/lp\. In the ultrarelativistic limit the uncertainty in the energy of the particle is AE ~ cAp, which is around 
one Planck energy 

£pi:=J^~2.0-10 9 J. (1.2) 

V Cr 

In other words, we have enclosed one Planck energy inside a cube whose diameter is one Planck length. That amount of 
energy, however, is enough to shrink the spacetime region bounded by the cube into a black hole with a Schwarzschild 
radius equal to around one Planck length. So it seems possible that when probing the structure of spacetime at the 
Planck length scale one encounters with Planck size black holes. 

At the present stage of research a construction of a precise mathematical model of spacetime as a whole made of 
tiny black holes is out of reach. However, it is possible to test the idea of spacetime as a foam of microscopic black 
holes in the context of an important specific problem of quantum gravity. The problem in question is the microscopic 
origin of black hole entropy. There is a general agreement between the researchers working in the field of quantum 
gravity that black hole has an entropy which is one quarter of its event horizon area or, in SI units, 
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This result is sometimes known as the Bekenstein-Hawking entropy law. In addition to black holes, it is valid for 
cosmological, de Sitter and Rindler horizons as well, and it seems to imply that in addition to the three classical, 
there is an enormous amount of quantum-mechanical degrees of freedom in black holes. It is reasonable to expect 
that these additional degrees of freedom lie at the horizon of the black hole. During the recent ten years or so several 
attempts have been made by experts of string theory and loop quantum gravity to identify the quantum-mechanical 
degrees of freedom of black holes, and to provide a microscopic explanation for the black hole entropy Q. 

To test the viability of an idea of spacetime as a foam of Planck size black holes we postulate, in this paper, a 
specific spacetime foam model of the event horizon of a Schwarzschild black hole |3j . It follows from our model that the 
entropy of the Schwarzschild black hole is proportional to its horizon area, and we also present arguments suggesting 
that the constant of proportionality must be, in natural units, equal to one quarter. In other words, it seems possible 
to obtain Eq. I|1.3fl from our model. 

This paper is organized as follows: In section 2 we introduce our model, especially the postulates on which it is 
based. In short, the event horizon of a macroscopic Schwarzschild black hole is assumed to consist of Planck size, 
independent black holes. Each of the microscopic holes on the horizon is assumed to obey a sort of "Schrodinger 
equation" of black holes, which was published some years ago. The postulates of our model connect the quantum 
states of individual microscopic holes on the horizon with the area eigenvalues of the event horizon of the macroscopic 
Schwarzschild black hole. 

In section 3 we proceed to calculate the entropy of the Schwarzschild black hole from the postulates of our model. 
An essential ingredient of our calculation of black hole entropy is our decision to apply classical statistics to our 
model. This decision of ours is motivated by Hawking's result that, at least in the semiclassical limit, black hole 
radiation spectrum is the continuous black body spectrum. If one attempted to apply any sort of quantum statistics 
to our model, a discrete radiation spectrum radically different from the one predicted by Hawking would follow even 
for macroscopic black holes. Continuous spectrum, however, is regained if classical, instead of quantum, statistics 
is applied to our model. We find that our model supplemented with classical statistics implies that in a very low 
temperature the entropy of the Schwarzschild black hole is, up to an additive constant, proportional to its horizon 
area. In section 4 it is claimed, on grounds of certain geometrical arguments, that it is natural to take the constant of 
proportionality to be equal to one quarter. This result reproduces the Bekenstein-Hawking entropy law of Eq. 1)1. 3|) . 

Section 5 contains a critical analysis of our model. We list the reasonable objections against our model we have 
managed to find, and our answers to these objections. 



The starting point of our model is Bekenstein's proposal of the year 1974, which states that the eigenvalues of the 
area of the event horizon of a black hole are of the form: 



where n is an integer and 7 is a pure number of order one. In other words, Bekenstein's proposal states that the event 
horizon area of a black hole has an equal spacing in its spectrum p| . 

Bekenstein's proposal has been revived by several authors on various grounds |5|. One way to obtain Eq. i|2.1|) 
for Schwarzschild black holes is to consider the following eigenvalue equation for the Schwarzschild mass M of the 
Schwarzschild black hole (unless otherwise stated, we shall always use the natural units where h = c = G = ks = I)'- 



This equation which, in a certain very restricted sense, may be viewed as a sort of "Schrodinger equation" of the black 
hole from the point of view of a distant observer at rest with respect to the hole, was obtained in Ref. [|J from the 
first principles by means of an analysis which was based on Kuchaf 's investigation of the Hamiltonian dynamics 
of Schwarzschild spacetimes. In Eq. I|2.2|) . a is the throat radius of the Einstein-Rosen wormhole in a foliation where 
the time coordinate at the throat is the proper time of an observer in a free fall through the bifurcation two-sphere. 
■0(a) is the wave function of the hole, and the real number s determines the inner product between the black hole 
states such that 



II. THE MODEL 
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In other words, the Hilbert space is chosen to be the space L 2 (R + , a s da) 

It was shown in Ref. [|| that the spectrum of M given by Eq. I|2.2[) is always discrete, bounded from below, and 
can be made positive by means of an appropriate choice of the self-adjoint extension. Moreover, it was shown that if 
s = 2, then for large n the WKB eigenvalues of the mass M are of the form 

M„ WKB = y^n+T. (2.4) 

Surprisingly, this WKB approximation for large n provides an excellent approximation for the mass eigenvalues even 
when n is small (i.e. of order one): At the ground state where n — 0, the result given by Eq. (|2.4[1 differs around 
one percent from the real ground state mass eigenvalue, and the difference very rapidly goes to zero when n increases 
(see Appendix A). So we see that because the event horizon area of a Schwarzschild black hole with mass M is: 

A = 16ttM 2 , (2.5) 

Eq. (|2.2|1 implies, as an excellent approximation, the following spectrum for the horizon area: 

A n = 327r(n+i). (2.6) 

In other words, we recover Bekenstein's proposal of Eq. I|2.1|) with 7 = 32-7T. 
We are now prepared to state the postulates of our model. Our postulates are: 

1. The event horizon of a macroscopic Schwarzschild black hole consists of independent, Planck size Schwarzschild 
black holes. 

2. Each hole on the horizon obeys Eq. (|2.2|l with s = 2. 

3. Holes in the ground state where n = do not contribute to the horizon area. 

4. When the hole is in the nth excited state, it contributes to the horizon an area which is proportional to n. 

5. The total area of the horizon is proportional to the sum of the areas contributed by the holes on the horizon. 

No doubt it is rather daring to assume that the microscopic holes on the horizon are independent of each other. 
However, at this stage of research that assumption is necessary if we want to make progress. Moreover, one may 
speculate with a possibility that when the density of black holes on the horizon is constant, the effects of all the rest 
of the holes to an individual hole on the horizon may cancel each other. Another bold assumption is contained in 
the Postulate 3. As such, the Postulate 3 views observable black holes as excitations of black holes in a ground state. 
When compared to the Postulates 1 and 3, the Postulates 2, 4 and 5 appear as rather natural. The Postulates 4 and 
5 imply that the horizon area is of the form 

A = a(n 1 +n 2 + ... + n N ), (2.7) 

where n\, n%, . . . , njsr are the quantum numbers associated with the area eigenstates of the holes on the horizon, and 
a is an unknown constant of proportionality. We shall consider the postulates of our model in more details in section 
5. 



III. ENTROPY 

When the concept of entropy of a black hole was introduced by Bekenstein and Hawking 30 years ago, all the 
derivations of the expression l|1.3fl for black hole entropy were performed by means of semiclassical arguments: When 
Hawking obtained an expression T = for black hole temperature yielding the expression (|1.3|l for black hole 

entropy, he treated spacetime classically and matter fields quantum-mechanically 8J. Later, in 1977, when Gibbons 
and Hawking calculated the black hole entropy by means of Euclidean path integral methods [t| , they simply used a 
semiclassical approximation to the path integral. In other words, the expression (|1.3|l for the black hole entropy is, to 
the greatest possible extent, a "semiclassical entropy" which we should obtain from a microscopic model of spacetime. 

At this point it is useful to recall what we actually mean by the very concept of entropy in quantum and classical 
physics: In quantum statistics the entropy of a system in a given macrostate is the natural logarithm of the number 
of microstates corresponding to that macrostate, whereas in classical statistics entropy is the natural logarithm of 
the phase space volume corresponding to the given macrostate ^(J- ^ may be shown that the entropy of quantum 
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statistics reduces to the entropy of classical statistics when the spectra of observables are not assumed to be discrete 
but continuous. Which statistics — classical or quantum — should we use for spacetime itself? 

The answer to this question depends on whether we consider the radiation spectrum of a black hole continuous or 
discrete. If the radiation spectrum is discrete, then the horizon area spectrum of the macroscopic hole, as well as the 
spectra of the microscopic holes on the horizon, is discrete and we must use quantum statistics. Consequently, if the 
radiation spectrum is continuous, the area spectrum is also continuous and we must use classical statistics. According 
to Hawking, the radiation spectrum is purely thermal, continuous black-body spectrum. Therefore the answer to 
our question is obvious: If we want to obtain for the black hole entropy an expression which is in agreement with 
Hawking's radiation law we must, of course, use classical statistics for the black hole itself. 

This conclusion has nothing to do with the question about whether we really believe black hole radiation spectrum to 
be continuous, as in Hawking's theory, or discrete as it is, for instance, according to Bekenstein's proposal. Curiously, 
it is possible to express an argument, based on Heisenberg's uncertainty principle, to the effect that when the effects 
of matter fields are so strong that they overshadow the quantum effects of s pac etime, then the discrete spectrum 
predicted by Bekenstein's proposal reduces to Hawking's black body spectrum Actually, this is what we assume 
here. More precisely, we assume that when the effects of matter fields are absent, the black hole area spectrum 
is discrete and follows Bekenstein's proposal but when the effects of matter fields are strong enough, the spectrum 
becomes continuous, and classical statistics may be applied to the black hole itself. 

So, let us calculate the classical entropy, or the natural logarithm of the classical phase space volume, of our model. 
Eq. (|2.2[) may be deduced from the classical Hamiltonian 

H=±-y + a*) 1 (3.1) 

of the Schwarzschild black hole (For the details see Ref. In this equation p is the canonical momentum conjugate 
to a. The numerical value of H is the Schwarzschild mass M of the hole. Therefore we find, using Eq. 12. 4|) . that for 
a single black hole with "mass" Mi and canonical coordinates a\ and p\ on the horizon, 



(p^+af) 2 = =2rn + l. (3.2) 



However, since we are now considering classical statistics, n± no more is an integer but it may be any non-negative 
real number. Moreover, it is no more possible to associate with Mi any sensible physical meaning as the "mass" of a 
hole on the horizon (A more detailed investigation of this issue is performed in section 5) . M\ simply is a parameter 
with the property that M\ is proportional to the area contributed by a single microscopic hole to the total horizon 
area. 

Eq. I|3.2fl now implies that 

^^+«?) 2 + - + 4^ +a ^ 2 = f + i (3 - 3) 

It can be shown (See Appendix B) that for fixed A and N this is a compact (2N — l)-dimensional hypersurface in a 
2iV-dimensional phase space. Its (2N — l)-volume, in the units of (2ttK) n ~ 1 / 2 , is of the form 

AT-l/2 

il = C(N) \N+—\ , (3.4) 




where C(N) is a function which depends on N only. For fixed N, C(N) is finite. 

At this point we should pause for a moment. In section 2 we threw away with great pomp and ceremony the vacuum 
contribution to the black hole horizon area by stating that holes in the ground state do not contribute to the horizon 
area. Eq. I|3.3|) . however, involves the term y which is equal to the vacuum contribution to the horizon area. Are we 
now quietly returning back the vacuum contribution to our calculations? 

The answer to this question is emphatetically no. The reason why the term ^ suddenly appears in Eq. (|3.3|) is 
that in Eq. I|3.3|) we calculate the phase space volume of our model, and the role of the term ^ is to take into account 
the vacuum contribution to the phase space volume. In other words, although we assume that ground states do not 
contribute to the horizon area, they nevertheless contribute to the phase space volume. Actually, as we shall see in a 
moment, in very low temperatures most of the phase space volume is really occupied by the unobservable vacuum, or 
black holes in the ground state. 

The entropy of the Schwarzschild black hole is the natural logarithm of the phase space volume Q corresponding 
to a fixed horizon area A: 



5 = Inn = lnC(A0+ IN - - In \N + — . (3.5) 
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Using Eq. (|2.5(l and keeping N as a constant we may obtain the temperature T of the Schwarzschild black hole: 

1 ( 0S\ AT-i 64ttM 




or, since 2V is assumed to be very large, 



It is interesting that if N, the number of microscopic holes on the horizon, is sufficiently small when compared to 
the mass M of the macroscopic hole then the temperature increases, instead of decreasing as in Hawking's theory of 
black hole radiation, as a function of M . 

At this point we introduce a new assumption to our theory. We assume that most microscopic holes on the horizon 
are in the ground state, where n = 0. More precisely, we assume that the average 

m + . . . + n N A 

n := = (3.8) 

N Na y ' 

of the quantum numbers rz-i, • • • , n^v has the property 

n<l. (3.9) 
Since A is proportional to M 2 , we find that this condition may be written as 

N M y ' 

When the assumption 1|3.9[) is employed, the second term on the right hand side of Eq. I|3.7|) vanishes, and the 
temperature becomes to: 

The assumption l|3.9[l is sensible when we consider the thermodynamics of a macroscopic black hole in a very low 
temperature: In a very low temperature the constituents of any system tend to be as close to the ground state as 
possible. The same is true also for our model: We see that the minimum temperature is achieved when the second 
term on the right hand side of Eq. 13.7fl . and therefore n, becomes as small as possible. 

From Eq. 1 If) it now follows that, up to an additive constant which depends only on the number N of the 
microscopic holes on the horizon, the entropy of the Schwarzschild black hole is, in a very low temperature, 

9 A 

S=—. (3.12) 

a 

In other words we have obtained, up to an undetermined constant of proportionality, the Bekenstein-Hawking entropy 
law of Eq. Ijl.3|l . Here we have considered Schwarzschild horizon only but our consideration may easily be generalized 
for other horizons as well. 



IV. CONSTANT OF PROPORTIONALITY 



It only remains to fix the constant a in Eq. I|3.12[l . The constant a was defined in Eq. (|2.7|) . From that definition 
it follows that if we know how the total area of the horizon depends on the areas of the holes on the horizon, we 
may calculate a. A comparison of Eqs. I|1.3|) and (|3.12|) yields the result that, in natural units, we should have 
a = 8, and our task is to obtain this value by means of geometrical arguments. Of course, nobody really knows what 
the spacetime geometry at the Planck length scale really looks like, and therefore one should take such arguments 
with a pinch of salt. However, it is possible to express arguments which, when not understood as a real description 
of spacetime structure at the Planck length scale but rather as a sort of statistical average of the behaviour of the 
microscopic holes on the horizon, may at least provide a heuristic aid of thought. 

To begin with, consider Fig. ^ In that figure we have drawn microscopic black holes on the horizon such that, at 
this stage, they do not overlap each other. A distant observer does not see a hole on the horizon as a sphere but as 
a circle with radius equal to its Schwarzschild radius Rs- The area of this circle is 7ri?|., which is one quarter of the 
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FIG. 1: In our model spacetime is constructed from Planck size black holes. A faraway observer sees each particular hole as 
a circle with a Schwarzschild radius Rs- When a hole on the horizon exactly fits inside a square, the area of the square is — 
times the area of the corresponding circle. 



horizon area of the corresponding microscopic hole. Eq. (|2.tj|) . together with the postulate 3, therefore implies that 
the area covered by a microscopic hole in the nth excited state is 87m. 

When considering Fig. Q however, one finds that the holes on the horizon do not cover the whole horizon but 
certain voids remain. A better way to cover the horizon is to use squares drawn around the holes such that each 
hole exactly fits inside the corresponding square. In that case the area of each square is ^ times the area of the 
correponding circle. Summing the areas of the squares we find that the horizon area is 



from which it follows that a = 32. In other words, a is too large with a factor of 4. 

A remedy for this problem may be found if we assume that the holes on the horizon overlap each other. Recall that 
we are obtaining an expression for the maximum entropy of a macroscopic Schwarzschild black hole. The maximum 
entropy may be gained if there is a maximum amount of microscopic holes on the horizon. That is because if two 
black holes on the horizon come together to form a single black hole, the area of the resulting black hole is larger 
than is the sum of the areas of the original black holes. In other words, a larger amount of the horizon is covered 
by a smaller number of degrees of freedom. Since the maximum entropy is achieved when the number of degrees of 
freedom is as large as possible, it follows that when two black holes on the horizon overlap in the state of maximum 
entropy, they still must not form a single black hole. 

The maximum overlapping between two holes with equal sizes is achieved when the center points of the holes lie 
on each other's event horizons (See Fig. |5J). In that case the holes are not yet "swallowed" by each other, and the 
distance between the center points of the holes is their common Schwarzschild radius Rs- As one can see from Fig. 
121 however, covering the horizon with black holes overlapping each other in this manner is equivalent to covering the 
horizon with holes with radius ^Rs and not overlapping each other. Therefore one finds that the horizon area now 
takes the form: 



A = 32(n 1 + n 2 + ... + n N ), 



(4.1) 



A = 8(ni + n 2 + . . . + tin), 



(4.2) 



which implies the desired result 



(4.3) 



When this result is substituted in Eq. (|3.12(l . we get 




(4.4) 



or, in SI units, 



S= - 



lk R c 3 



A. 



(4.5) 



4 hG 



7 




I / 

FIG. 2: Two identical black holes on the horizon overlapping each other such that the center points of the holes lie on each 
other's event horizons. Covering the horizon with black holes overlapping each other in this manner is equivalent to covering 
the horizon with non-overlapping holes having radiuses which are one halves of the radiuses of the original holes. 



In other words, we have exactly recovered the Bekenstein-Hawking entropy of Eq. 

However, it must again be strongly emphasized that the geometrical arguments given here should not be considered 
as a real description of the structure of spacetime at the Planck length scale. Rather, they should be viewed as mere 
heuristic aids of thought. The microscopic structure of spacetime is certainly very much more complicated than is 
the simple picture given in this section. Our geometrical arguments, however, have a certain statistical content which 
may be expressed as follows: The average number density of the microscopic black holes on the horizon is such that, if 
the holes were understood as classical objects, their center points would lie, in average, on each other's event horizons. 
Whether one accepts this argument or not is, at the present stage of research, a matter of taste but nevertheless it 
produces the desired value for the constant a. 

V. SOME OBJECTIONS 

As it may have become clear to the reader, our model and the postulates on which it is based, contain several ideas 
which are physically new. Always when expressing a new physical idea it is of vital importance to try to prove oneself 
wrong, and this is what we shall attempt to do in this section. More precisely, we list all the reasonable objections 
against our model we have managed to find, and our answers to these objections. Our objections are: 

1. It is very strange to think a black hole horizon of being made of Planck size black holes. If black holes are made 
of Planck size black holes, then what are the Planck size black holes made of? 

Answer: In our model we have assumed the Planck size black holes described by Eq. (|2.2|) to be fundamental 
constituents, or building blocks, of spacetime. So they are assumed to have no intrinsic structure. The idea that 
spacetime might be made of Planck size black holes in one way or another is supported, among other things, 
by the arguments based on Heisenberg's uncertainty principle expressed in the Introduction. So far we have no 
precise model of the whole spacetime being made of Planck size black holes, and a construction of such a model 
should be viewed as a challenge for the future. We only have a model of a Schwarzschild horizon and, as one 
can see, our model reproduces the Bekenstein-Hawking entropy law. 

2. Although the paper intends to provide a model of a Schwarzschild horizon, an assumption that a horizon of 
spacetime is considered, is nowhere explicitly used in the derivation of the Bekenstein-Hawking entropy law. 
Actually it seems that arguments similar to those employed in this paper may be used to imply that every 
spacelike two-surface of spacetime has an entropy which is proportional to its area, no matter whether that two- 
surface is a horizon or not. Is this not in a contradiction with the known properties of gravitational entropy? 

Answer: It is true that an assumption that a horizon of spacetime is under scrutiny is nowhere used in this paper. 
The primary reason for our choice to construct a model of a Schwarzschild horizon is that with a Schwarzschild 
horizon one may associate a well-defined concept of energy which, in turn, may be used when calculating the 
temperature of the horizon. It is also true that one could argue by means of reasoning similar to that used 
in this paper that entropy is associated not only with spacetime horizons but with every spacelike two-surface 
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of spacetime. More precisely, it really seems that every finite spacelike two-surface, no matter whether that 
two-surface is a horizon or not, possesses an entropy which, in natural units, is one quarter of its area. 

Although astonishing, this conclusion, however, does not contradict with the known physics. Actually, there 
seems to be very good grounds to believe that every piecewise smooth, spacelike two-surface of spacetime indeed 
carries a certain amount of entropy although that entropy produces observable physical effects only for the 
observers having that two-surface as a horizon. Provided that one accepts the view held at least implicitly by 
many authors that horizon entropy is due to the microscpic degrees of freedom of spacetime at the horizon 
this result is a direct consequence of the well-established result that any finite part of a Rindler horizon has an 
entropy which, in natural units, is one quarter of the area of that part. The line of reasoning producing this 
conclusion may be expressed as follows: 

(a) Each finite part of a Rindler horizon possesses an entropy which is one quarter of its area. 

(b) Each finite spacelike two-plane is, from the point of view of an appropriate observer, a part of a Rindler 
horizon. 

(c) Therefore, each finite spacelike two-plane possesses an entropy which is one quarter of its area. 

(d) Each piecewise smooth spacelike two surface is a union of infinitesimal spacelike two-planes, each having 
an entropy equal to one quarter of its area. 

(e) Therefore, each piecewise smooth spacelike two-surface possesses an entropy which is one quarter of its 
area. 

In other words, an explicit assumption that one looks at the horizon is not needed when one attempts to obtain 
the Bekenstein-Hawking entropy law: There are good grounds to believe that the Bekenstein-Hawking entropy 
law holds for any piecewise smooth, spacelike two-surface. We hope to be able to consider this subject in more 
details in the forthcoming papers. 

3. How can the Planck size black holes on the horizon be independent of each other? Certainly there should be an 
extremely strong interaction between them. 

Answer: This is indeed what one might expect on classical grounds. However, we are now considering spacetime 
at very small length scales where quantum gravity reigns, and all bets are on. It is possible that the interactions 
between black holes effectively cancel each other in a bit similar way as the gravitational effects cancel each 
other from the point of view of an observer at the center of the Earth. This is an open question. 

In addition to these rather general remarks it is possible to express a much more serious argument supporting 
the claim that the microscopic black holes might indeed be independent of each other. This argument is based on 
Jacobson's exremely important discovery that Einstein's equations may be viewed as thermodynamical equations 
of state |l2j ■ More precisely, Jacobson showed that if one assumes that a local Rindler horizon always possesses 
an entropy which is one quarter of its area, then Einstein's equations follow from the first law of thermodynamics. 
So one is inclined to think that Einstein's equations, with their prediction that macroscopic black holes attract 
each other, are probably not very fundamental at all but they are merely a consequence from the statistics of 
the fundamental constituents of spacetime. There is not necessarily any interaction between these fundamental 
constituents at the microscopic level, but their statistical properties imply at the macroscopic level properties 
which might lead one to conclude an existence of a sort of effective interaction. As a familiar example one 
may think of a classical ideal gas: Classical ideal gas has a certain pressure which might lead one to conclude 
(incorrectly) that between the particles of the gas there is a repulsive interaction which prevents one to compress 
the gas. However, no such repulsive interaction really exists, and the pressure is simply a consequence from 
the statistical mechanics of the gas. Perhaps something similar happens with microscopic black holes: At the 
microscopic level there is no interaction between the holes but their statistical properties imply in a certain 
limit, among other things, Einstein's equations with all their predictions. Actually this may well be the case, 
because it seems to us that our analysis could be generalized to show that not only black holes but also a 
local Rindler horizon possesses an entropy which is one quarter of its area, and, as Jacobson has showed, a 
generalization of the Bekenstein-Hawking entropy law for local Rindler horizons implies Einstein's equations. 
So it is possible that the gravitational interaction between macroscopic black holes may be a direct consequence 
from the independence of microscopic black holes. 

4. Is it not absurd that black holes in the ground state where n = do not contribute to the horizon area? 

Answer: This is probably the most serious objection against our model. However, the idea that vacuum states 
do not produce measureable effects (except in very special circumstances) is not very unfamiliar in physics: In 
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quantum field theories in flat spacetime the field is decomposed into its Fourier components, and the possible 
energies of each component are of the form: 



k.n 



(n+i)/*^, (5.1) 



where n = 0, 1, 2, • • • and is the angular frequency of the component corresponding to the wave vector k. 
The vacuum energy where n — 0, is never observed but only the energies which are above the vacuum energy 
may be measured. In other words, everything one can observe are excitations of the vacuum but the vacuum 
itself cannot be observed (except by means of the Casimir effect). It is not entirely impossible that this feature 
of nature might pertain also to quantum gravity. In quantum gravity, however, energy is replaced by area: As 
one can see from Eq. (|2.(j[) . the possible areas of the microscopic holes on the horizon are exactly of the same 
form as are the energies of the Fourier components of a quantized field. One is therefore tempted to draw an 
analogy between the quantization of energy in ordinary quantum field theories, and the quantization of area in 
our model: In the same way as everything one observes in ordinary quantum field theories arc excitations of the 
vacuum state of energy, in our model everything one observes are excitations of the vacuum state of area. This 
analogy may be viewed as the main justification for the Postulate 3. 

5. It follows from the Postulate 5 that the squares M 2 of the masses of the individual holes on the horizon sum 
up, essentially, to the square M 2 of the total mass of the hole. Is this not wrong because certainly the masses 
themselves should sum up to the total mass of the hole? 

Answer: The concept of mass is very problematic already in classical general relativity. Actually, the concept 
of mass may be properly defined only when spacetime is asymptotically flat, and in that case the relevant mass 
concept is the ADM mass, which measures the mass-energy included by the whole spacetime. When we try to 
define the mass-energy included by a specific part of curved spacetime, however, we run into grave difficulties 
because the energy-momentum tensor of the gravitational field, in general, is not well-defined. So it is not clear 
what we are talking about when we say that the mass of a certain object in curved spacetime is such-and-such, 
unless that is the only object in an asymptotically flat spacetime, and we may define its ADM mass. 

In particular, the concept of mass of an object at the event horizon of a black hole is very problematic. To see 
how problematic it really is, consider the four-momentum of a particle having a rest mass m in flat spacetime: 

pf* = mf, (5.2) 

where the dot means proper time derivative. When spacetime is static, i.e. the metric tensor does not depend 
on the time parameter x , the quantity 

Po = mg 0fi x^ (5.3) 

is conserved along geodesies and also along curves having timelikc Killing vectors as their tangents, and it offers 
the best attainable definition for the concept of energy of a particle in curved spacetime. Since the Schwarzschild 
metric may be written as 



ds 2 - (l - — )dt 2 - - r 2 d9 2 - r 2 sin 2 9 def , (5.4) 

^ r 

we find that 

p =m(l-—)i (5.5) 



r / 

For a particle at rest with respect to the Schwarzschild coordinates we have 



dT 

and therefore 

2M\V2 



. dt / 2M\-V2 

t = — = {\ , (5.6) 



pa = mil , (5.7) 



r i 



Hence we see that when an object lies at the horizon where r = 2M, its energy, according to above mentioned 
definition, is zero. Of course it is impossible to keep a particle still at the horizon but nevertheless our consider- 
ations provide an example of the problems of the mass-energy concept at the event horizon of a black hole, even 
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in the classical level. In the quantum level one may expect even deeper problems. Therefore, if we consider the 
event horizon of a black hole as a system made of tiny black holes, the mass parameters M, of the holes do not 
represent mass in any ordinary sense. Consequently, the sum of the parameters Mj is not the Schwarzschild mass 
M of the macroscopic hole, either. However, each hole on the horizon covers a region with an area proportional 
to Mf on the horizon. Therefore it is natural to postulate that the sum of the quantities Mf equals, up to a 
constant of proportionality, to the square M 2 of the Schwarzschild mass of the hole. 

6. The entropy of the hole has been calculated by using classical statistics right from the beginning. Would it not 
have been more appropriate to perform the calculations by using quantum statistics first, and then take the 
classical limit? 

Answer: This question has been considered in details already in section 3. Still, it is appropriate to repeat here 
the essential points: It follows from Eq. (|2.7(l which, in turn, is a direct consequence from our postulates, that 
the radiation spectrum of a black hole is discrete. However, according to Hawking's radiation law the black 
hole radiation spectrum is a continuous black body spectrum. It is impossible to obtain Hawking's continuous 
spectrum from our model unless one assumes right from the beginning that the quantum numbers nj., »i2, • • • , njy 
may take any values, and classical statistics applies. It follows from Eq. Ij2.7|l that the possible frequences of 
the quanta of radiation are, for a = 8, integer multiples of the fundamental frequency 

Even for macroscopic black holes (i.e. when M is of the order of ten solar masses) this quantity is fairly big: It is 
of the order of 0.3 kHz which is about the same as is the resolving power of an ordinary portable radio receiver! 
So if one wants Hawking's continuous spectrum out from the discrete spectrum predicted by Eq. Q2.7|). the 
stationary quantum states of the black holes must become very much mixed with each other such that classical 
statistics, in effect, may be applied. A mechanism for this kind of mixing through the interaction between the 
hole and matter fields was proposed in Ref. fill ]. 

7. The expression of Eq. I|3.5|l for the black hole entropy is not exactly one quarter of the horizon area, even for big 
JV, but it contains an additive constant which depends only on TV. The presence of that term should produce 
physical effects different from those obtainable from the Bekenstein-Hawking entropy law. 

Answer: Actually, the quantity N represents the particle number in our model, the "particles" being now the 
microscopic black holes on the horizon. In all thermodynamical systems the physically observable thermody- 
namical quantities are related to those partial derivates of entropy where the particle number N is kept constant 
during the differentiation. For instance, the temperature T of a system is the inverse of the partial derivative 
of entropy S with respect to energy E such that the particle number N and the volume V of that system are 
kept constant: 



1 

f 




(5.9) 

V,N 



Moreover, the pressure p of a system is its temperature times the partial derivative of entropy with respect to 
the volume V such that energy E and particle number N are kept constant: 



p--=A d £\ ■ (5.10) 



\dV 



E.N 



The only quantity which in any thermodynamical system is related to the partial derivative of entropy with 
respect to the particle number TV is the chemical potential 



M:=-r ™ ( 5 - n ) 



E,V 



of the system. As far as one is not interested in the chemical potential of the Schwarzschild black hole, the 
physical predictions obtainable from our expression for black hole entropy are identical to those obtainable 
from the Bekenstein-Hawking entropy law. As it comes to the chemical potential of the Schwarzschild black 
hole, in turn, nobody has the slightest idea about what that might be for the very simple reason that nobody 
knows what the "particles" constituting a black hole would be. So we see that our expression for the black 
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hole entropy does not contradict with the existing knowledge ol black hole thermodynamics. In general, it 
seems possible, at least in principle, to use different chemical potentials to distinguish different derivations of 
the Bekenstein-Hawking entropy law. More precisely, even if a certain derivation of the Bekenstein-Hawking 
entropy law yielded an expression for the temperature of the horizon consistent with that law, the number of 
degenerate states corresponding to the same horizon area A is not necessarily e^ A but an expression for the 
number of degenerate states may contain a model-dependent prefactor, which is a function of the chemical 
potential of the system. 

VI. CONCLUDING REMARKS 

In this paper we have considered a spacetime foam model of the Schwarzschild horizon where the horizon of 
a macroscopic Schwarzschild black hole consists of Planck size Schwarzschild black holes. Using this model we 
found that the entropy of a macroscopic Schwarzschild black hole is, up to an additive constant, proportional to 
its horizon area. It seems that our derivation of this result is valid for any horizon but here we have restricted 
our consideration to Schwarzschild horizons only. Certain heuristic arguments may be employed to imply that the 
constant of proportionality is, in natural units, equal to one quarter. 

The key points in our derivation were an assumption that when a microscopic black hole on the horizon is in a 
ground state, it does not contribute to the horizon area, and our decision to apply classical statistics to our spacetime 
foam model. An assumption that a hole in a ground state, even if its Schwarzschild mass were non-zero, does not 
contribute to the horizon area, may be motivated by an analogy with ordinary quantum field theories: In ordinary 
quantum field theories the vacuum state of energy is not observed, whereas in our model the vacuum state of energy is 
replaced by the vacuum state of area. As it comes to our decision to apply classical, instead of quantum statistics to 
our model, that decision may be considered justified on grounds of the fact that according to Hawking the radiation 
spectrum of a black hole is a continous black body spectrum. If we want to obtain a continuous spectrum for that 
radiation we must assume that the spectra of observables of a macroscopic black hole are continuous, and therefore 
we must use classical statistics for the black hole itself. 

In our model we assumed at first that the tiny holes on the horizon obey Bekenstein's proposal. In other words, we 
assumed that when the effects of matter fields are neglected, the horizon area spectra of the holes on the horizon have 
an equal spacing. That spacing was obtained from Eq. 1)2.2)1 . the "Schrodinger equation" of the Schwarzschild black 
hole. When the effects of matter fields are assumed to be strong enough, however, the discrete area spectrum washes 
out into continuum, and Hawking's continuous black body spectrum is recovered. In other words, we considered the 
implications of Eq. 12.21 in a limit where the uncertainties in the area eigenvalues are of the same order of magnitude 
as is the spacing between the area eigenvalues of nearby states. The fact that we used Eq. l|2.2[l in our derivation 
of an expression for the black hole entropy may be viewed as an argument supporting the physical validity of that 
equation as well as of Bekenstein's proposal. In this paper we do not express opinions about whether the radiation 
spectrum really is continuous or discrete. We only showed that if the radiation spectrum is assumed to be continuous, 
then the Bekenstein-Hawking entropy law follows from our model. 

Taken as a whole, our model may be viewed as an attempt to understand the microscopic structure of spacetime. 
No doubt the picture provided by our model is very different from those provided by, for instance, string theory and 
loop quantum gravity. An andvantage of our model, however, is that it takes seriously the possibility suggested, 
among other things, by Heisenberg's uncertainty principle, that spacetime at the Planck length scale might consist 
of Planck size black holes. So far we have no precise model of spacetime itself but only of its Schwarzschild horizon. 
However, our results with the derivation of the Bekenstein-Hawking entropy law are encouraging, and it will be very 
interesting to see whether the ideas presented in this paper may be worked out into a precise mathematical model of 
the microscopic structure of spacetime. 
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APPENDIX A: MASS EIGENVALUES 



During the preparation of this paper we have studied the mass eigenvalues of Eq. (|2.2() . when s = 2, for small n in 
two different ways. The first method we have used is perturbation theory. The second method is to solve Eq. I|2.2|) 
numerically. Both of these methods, of which the second one is more reliable, seem to imply that Eq. I|2.4|l provides 
an excellent approximation for the mass eigenvalues even for small n. 



Perturbation Theoretic Approach 



As a starting point we consider Eq. 12.21) when s = 2: 




We write the Hamiltonian of Eq. ijAlfl in the form 



H 



1 d 2 



1 d 1 



2a da 2 2a 2 da 2 



%p{a) = Mtp{a) 



a=:Hn+H\ 



where 



and 



Hn ;= - — 



1 d 2 1 



2a da 2 2° 



H> :=- — -. 

2a 2 da 



If we neglect the term H' in Eq. (|A1|) and denote 



u := a - M, 



we see that the resulting differential equation takes the form: 



(Al) 



(A2) 



(A3) 



(A4) 



(A5) 



(A6) 



This equation is similar to that of a one-dimensional harmonic oscillator, and therefore members of one set of solutions 
to Eq. (|A6|I are of the form 



V4 0) (a) =N n H n {a-M^)e 



■i(a-M^)) 2 



(A7) 



where H n denotes Hermite polynomial of order n, N n is a normalization constant, and the eigenvalues Mn are 

M<°> = V2n+ 1. (A8) 



One can see that the eigenvalues of Eq. (| A8|) are the same as the WKB eigenvalues in Eq. (|2.4|l . 

Now we consider the term H' as a small perturbation of the unperturbed Hamiltonian Ho. Here we must use 
somehow "non-standard" methods for calculating the first order perturbation because the Hamiltonian Ho is not a 
hermitian operator with respect to the inner product 1)2.3(1 . Usually, when using perturbation theory, one writes the 
Hamiltonian in the form H = Ho + AiJ', where the perturbation parameter A e [0, 1]. During the calculation of first 

order perturbation one must operate with Hq to the bra- vectors (ipr? \ ■ The problem is that if Hq is not a hermitian 
operator, it cannot be transported from the right to the left hand side inside the Dirac brackets, and the first order 
perturbation cannot be calculated in a standard way. However, if we assume that our perturbation expansion is viable 
when A = 1, we can set right from the beginning that A = 1. In this case the operator Hq + XH is indeed a hermitian 
operator and one can operate with this operator to the bra-vectors in the usual manner. Using this trick, one can 
calculate the first order perturbation. Let us begin with the perturbation expansions 



|$„) = |^ 0) >+A|^)+A 2 |^ 2 



(A9) 
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M„ 



M(°) + AM« + A 2 M< 2 ) 



(A10) 



where the vectors |$ n ) are eigenstates of the Hamiltonian H and the numbers M n are the corresponding eigenvalues. 
Symbols \ipn) and M„ ' denote corrections of order k to the unperturbed eigenstate and to the eigenvalue 



Af! 0) . When these expansions are inserted into the equation (ipn '\H\<& n ) = (ipn |M n |$„), we get: 
(i>W\(H + XH')\^) + X(^\(H + XH'MP) + A 2 (ip^\(H + XH')\^) + ...= 

(^|^)+A<4 0) I^ 1) )+A 2 (^ )|^ 2 ) 



(0) 



'.<°>l 



(All) 



Mi°> + AM« + A 2 Mi 2 ) 



((i?o + XH')^\^) + \((H + XH'WPWP) + X 2 ((H + XH')^\^) 



+ AM« + A 2 M r (2) 



1 + X(^\^)+X\^\^) 



Hence the first order perturbation is of the standard form 

M« = <4 0) |£'|4 0) > 



(A12) 



Let us now calculate the first order perturbation explicitly by using the definition (|2.3|l of the inner product. We 



get: 



m« = mHW) = / ^ or ( a )(-^i^ o) («))« 



da, 



(A13) 



= ^ H n {a M(°)) e -MW [<( a - M <°>) + (M(°) - a) ■ H n (a- M<°)) 

where H' n (a — M^) = 4-H n (a — M„), and the coefficients N n are determined by the requirement 

It turns out that the integral in Eq. I|A13(I can be written in a very simple form: 

Integrals in Eq. (|A14(I can be solved analytically by means of the error function erf (a:), but unfortunately a simple 
general formula for the coefficient N n seems to be somehow out of reach. However, one easily finds that 



(A14) 



(A15) 



TV 2 



* e -i + 3^(1 + erf (1)) 



and therefore, according to Eq. I)A15|) . we have 



A I, 



(i) 



| + | eV5F(l + erf (1)) 



-i -l 



0.0349228. 



(A16) 



(A17) 



First order perturbations in the cases where n = 1,2, 3, . . . can also be evaluated analytically but integration in 
Eq. l|A14fl becomes more complicated when n increases. Therefore, one might find some mathematical programs (like 
Mathematica) as very useful tools when one tries to calculate the coefficients N n analytically. In Table [I] we give the 
numerical values of first order perturbations in the cases where n = 0, 1, . . . , 10. 

However, the perturbation method has two weak points. Firstly, we cannot be sure that our perturbation expansion 
converges. We can only hope that first order perturbations provide a good approximation to the difference between 
the real and the WKB eigenvalues of the Schwarzschild mass M. Secondly, it is hard to find out to which self-adjoint 
extension our perturbation approach is connected (For details see Ref. @). For these reasons our results for first order 
perturbations can be thought only as suggestive and we need a more reliable method for the study of the eigenvalues 
of M. This method will be discussed in the next section. 
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n 


M< 1J 


M n 


m wkb 





0.0349228 


1.01 


1 


1 


0.0093795 


1.74 


1.7320 . 


2 


0.00513602 


2.24 


2.2360 . 


3 


0.00346 


2.65 


2.6457. 


4 


0.00257746 


3.01 


3 


5 


0.00203793 


3.32 


3.3166. 


6 


0.00167622 


3.61 


3.6055 . 


7 


0.001418 


3.88 


3.8729 . 


8 


0.00122503 


4.13 


4.1231 . 


9 


0.00107574 


4.36 


4.3588 . 


10 


0.000957045 


4.59 


4.5825 . 



TABLE I: Numerical values of M„ and M™ KB . The results are given with the best possible precision our numerical 

analysis may offer. One may observe that WKB eigenvalues for the mass M provide an excellent approximation to the exact 
mass eigenvalues even when n is small. 



2. Numerical Analysis 



The numerical analysis of the mass eigenvalues has been performed with Femlab, and the results can be found in 
Table [I] In our numerical analysis we solve eigenvalues for the differential equation 

(-i»-5? + 5^*)*w- J, *w (A18) 

with the boundary conditions ip(0) = i/j (oo) = (Of course, the boundary condition ?/>(oo) = must be put in as 
ip(a) = 0, where a is a "sufficiently" large number). It may be shown that the mass spectrum given by this equation is 
identical to the mass spectrum given by Eq. IjAlj) Q. One can also study which self-adjoint extension the numerical 
routine has chosen when calculating the mass eigenvalues: The eigenfunctions corresponding to M n in Table^seem to 
behave, within the limits of the precision of the numerical computing, like ip(%) ~ y/x when x — » 0. This observation 
fixes the self-adjoint extension for our eigenvalues. 

Finally, it should also be noted that, due to a limited accuracy of numerical computing, some of the eigenvalues 
M n might be around one percent smaller than the values in Table [i] This does not cause any problems and actually 
the mass eigenvalues are brought even closer to the WKB estimate of Eq. I|2.4|l . 



APPENDIX B: PHASE SPACE VOLUME 



In our model, phase space volume corresponding to a fixed horizon area A is determined by the condition (|3.3|) : 

Pt+aj) +... + ^-n-(p 2 N +a 2 N ) =- + -. (Bl) 



i\ 1 8a^ 



a 



This equation describes a closed, compact (2N — l)-dimensional hypersurface E in the 2iV-dimensional phase space 
spanned by the canonical coordinates a, and pi. To calculate its volume we choose coordinates Ai, • • ■ , Aat_i G [0, 1] 
and <pi,"- ,<Pn € [0, 2tt] on the hypersurface S such that 

ai = L(Xi + Ai cost^i), (B2) 

Pi = LXi sint/Ji, (B3) 

a 2 = L(X 2 + A 2 cos</j 2 ), (B4) 

p 2 = LX 2 sin <f 2, (B5) 

a,N — L(Xn + Xn cos tpjsr), (B6) 

Pn = LXn sin ippf. (B7) 
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Here L — y N + and A^v = Ajv(Ai, • • • , \n~i) can be calculated from Eq. IjBljl : If one substitutes the coordinates 
a,i and pi in terms of Ai, • • • , \n and (p±, ■■ ■ , ipN into Eq. (|Blfl . one gets 



X N = ^l-Xl-Xl-...-X%_ v (B8) 

Eq. (|B8|) represents a hypersphere in the iV-space spanned by the coordinates A^ . So it is natural to perform another 
coordinate transformation from the coordinates Ai to the generalized spherical coordinates 9\, ■ ■ ■ , On -2 £ [0, V\ and 
0G[O,f]: 

JV-2 



A x = cos0 Y[ sinflj, (B9) 

N-2 

A 2 = sin0 J| sin6»,;, (BIO) 



i=l 
N-3 



A 3 = cos 6»i sin^, (Bll) 

i=l 
N-4, 

A 4 = cos 6» 2 sin (9,, (B12) 



i=l 



Ajy = cos0jv-2- (B13) 
Generally, the volume of the any (smooth) A^-dimcnsional hypersurface 5* can be evaluated from the integral 



J V9d N x, (B14) 



where g is the determinant of the metric. In our case the position vector of a given point on the hypersurface E can 
formally be written as 

r = ai(cpi, ■ ■ ■ ,ipn,Qi, ■ ■ ■ ,6n-2,4>)zx +Pi{ l P\, • • • ,Vn,8i, • • ■ ,0jv-2,<£) e 2 + . . . (B15) 
+ajv(vij ' ' ' 7 Viv, #i, • • • j 0N-2, <t>) ejv-i + Piv(</>1, • • • , fN,&i, ••• , On -2, 4>) ejv, 

where e,'s (i = 1, • • • ,2N) are orthogonal unit vectors in our 27V-dimensional phase space. From Eqs. i|B2|) - (|B7l) it 
follows that all non-vanishing components of the metric on the hypersurface, defined by the relation 

9ij = bi ■ bj, (B16) 

where bi — J^- (i = 1, ■ ■ ■ , 2N — 1) is the tangent vector of the coordinate curve corresponding to the coordinate a; 1 , 
are proportional to L 2 , and otherwise contain only products of the sines and the cosines of the angles (j>, Qi and v?,. 
Therefore the phase space volume ft corresponding to fixed area A has, in the units of (2nh) N ^ 1 / 2 , the form: 

tt/2 tt/2 tt/2 2tt 2tt 



2JV-1 



L 



n = ( 27 r)^-V2 J d<t> J ddl ' ' ' J d0N ~ 2 J d ^'"J ^Wtf. (B17) 
00 

where g is the determinant of the metric whose components g^ are obtained from the components gij such that we 
simply divide each g^ by L 2 . 

The (2N — l)-dimensional integral in Eq. (|B17(I is, unfortunately, very hard to calculate. However, one can easily 
see that this integral converges: From Eqs. (|B2(l - i|B7|) and l|B9|l - (|B13|) it follows that g is everywhere finite (it contains 
only products of the sines and the cosines of the angles 0, Bi and cpi). Furthermore, the integration is performed over 
a bounded region because <pi £ [0, 2ir] and 9i, <f> £ [0, -|]. Therefore, the (2N — l)-dimensional integral in Eq. (|B17|) 
converges. Its value, however, depends only on N and we can denote 

tt/2 tt/2 tt/2 2tt 2tt 

jL_ J dct> J dOx--- J d6 N -2 J d Vi---J d VN ^=:C{N). (B18) 



(2tt) n - 1 / 2 



00 
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So we can write the phase space volume f2 as: 

/ v N-l/2 

n = C{N)(N+^\ , (B19) 

which is Eq. I|3.4|) . 
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